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NOMENCLATURE,  SYMBOLS,  AND  ASSORTED  MARKS 


-  rectangu]ar  Cartesian  coordinates 

-  general  curvilinear  coordinates 

-  general  coordinates  on  middle  surface  of  undeformed 
shell 

1 

-  coordinate  normal  to  middle  surface  of  undeformed  shell 

-  polar  parameters  on  middle  surface  of  undefomed  shell 

-  cartesian  parameters  on  middle  surface  of  undeformed 
shell 

t 

-  radius  of  revolution  of  middle  surface  of  undeformed 
paraboloid 

e 

-  angular  coordinate  of  middle  surface  of  undeformed 
paraboloid 

r 

-  slope  of  meridian  tangent  to  middle  surface  of 
paraboloid 

-  angle  which  meridianal  tangent  makes  with  tangent  plane 
to  apex  of  paraboloid  (see  figure  2.1.3) 

-  focsil  length  of  middle  surface  of  paraboloid 

-  radius  vector  to  point  on  middle  surface  of  undeformed 
shell 

f 

-  radius  vector  to  point  in  undeformed  shell 

vli 


>^0 

-  radius  vector  to  point  on  middle  surface  of  deformed 
she.l] 

R 

-  radius  vector  to  point  in  deformed  she].) 

-  covariant  base  vectors  of  middle  surface  of  undeformed 
shell 

-  contravariant  base  vectors  of  middle  surface  of  un- 
deforned  shell 

h 

-  unit  normal  to  middle  surface  of  undeformed  shell 

fi  >  U’ t-3 

-  covariant  base  vectors  of  undeformed  body 

f 

-  contravariant  base  vectors  of  undeformed  body 

f'mn 

-  covariant  metric  tensor  of  undeformed  body 

-  m(\ 

-  contravariant  metric  tensor  of  undeformed  body 

^mn 

-  covariant  metric  tensor  of  middle  surface  of  undeformed 
shell 

^  m  A 

OL 

-  contravai'iant  metric  tensor  of  micMle  surface  of 
undeformed  shell 

covariant  base  vectors  of  deformed  body 

G^,  63 

-  contravariant  base  vectors  of  deformed  body 

A  ^  ,  A 

-  covari.ant  base  vectors  of  deformed  middle  surface  of 
shell 

viii 


A\  a2 

-  contravariant  base  vectors  of  deformed  middle  surface 
of  shell 

N 

-  unit  normal  vector  to  deformed  middle  surface 

Rj  1R2 

-  principal  radii  of  curvature  of  middle  surface  of 
undeformed  shell 

^11  >  ^22 

-  second  fundamental  quadratic  covariant  tensor  of 
undeformed  middle  surface  of  shell 

bj,b2,  bl,  bf 

-  second  fundamental  quadratic  mixed  tensor  of  undeformed 
middle  surface  of  shell 

y  y  m  y  rn  n 

^mo  ’  '’n  » 

-  strain  tensors 

-  stress  tensors 

gOL/i 

-  stress  resultant  tensor  of  shell 

-  moment  resultant  tensor  of  shell 

i“ 

-  transverse  shear  resultant  tensor  of  shell 

No(/3 

-  pliysical  components  of  force  resultants  tensors  referred 
to  5  \  C  ^  coordinate  system  (units  of  force  per  unit 
length) 

^rt/3 

-  physical  components  of  moment  resultants  tensors 

coordinate  system  (units  of  force-length  per  unit  length) 

<^a 

-  physical  components  of  transverse  shear  tensor  co- 

orciinate  system  (units  of  force  per  unit  length) 

h  »  if’  ^5 

1  2.  3 

-  unit  base  vectors  associated  witli  Y  ^  Y  >  / 

ix 


-  covariant  permutation  surface  tensor 

-  contravariant  permutation  surface  tensor 

-  force- resultants  referred  to  r,  ^  coordinate  system 
(units  of  force  per  unit  length) 

-  moment-resultants  referred  to  coordinate  system 

(units  of  force-length  per  unit  length) 

1  Qe 

-  transverse  shear  resultant  referred  to  r,  ^  coordinate 
system  (units  of  force  per  unit  length) 

-  focal  length  of  parabola 

-  powers  of  f,  the  focal  length 

.  <r 

CL  LL 

>  ^<r 

-  displacement  tensors  of  middle  surface 

^O- 

-  rotation  tensors  of  middle  surface 

or 

-  displacement  of  middle  surface  along  ^  (units  of 

length) 

0 

-  displacement  of  middle  surface  along  tangent  to  meridian 
(units  of  length) 

0 

-  displacement  of  middle  surface  along  tangent  to  latitude 
(units  of  length) 

w 

-  strain  tensor  of  middle  surface 

^a/S 

-  strain-curvature  tensor  of  middle  surface 

X 


extensional  strain  along  meridian  (dimensionless) 


oxtensional  strain  along  latitude  (dimensionless) 


shear  strain  (dimension.! ess) 


extensional  strain  of  middle  surface  along  meridian 
(dimensionless) 


extensional  strain  of  middle  surface  along  latitude 
(dimensionless) 


shear  strain  of  middle  surface  (dimensionless) 


extensional  strain-curvature  of  middle  surface  along 
meridian  (units  of  (length)"^) 


extensional  strain-curvature  of  middle  surface  along 
latitude  (  units  (length)  ) 

shear  strain-curvature  of  middle  surface  (units  of 
(length)  ) 

thickness  of  "undeformed  shell  (units  of  length) 
Youngs  modulus  (units  of  force  per  unit  area) 
shear  modulus  (units  of  force  per  unit  area) 

Poissons  ratio  (dimensionless) 
weight-density  (units  of  force  per  unit  volume) 


V  ^ 

-  powers  of  V  ,  Poissons  ratio 

-  coordinate  normal  to  middle  surface  of  deformed  shell 

-  thickness  of  deformed  shell  (units  of  length) 

1\1^ 

-  force-resultant  vectors  (units  of  force  per  unit  length) 

fri’ , 

-  moment- resultant  vectors  (units  of  force-length  per 
unit  length) 

-  tensor  components  of  body  force  vector 

P 

-  body  force  vector  (units  of  force  per  unit  area) 

’  f  ^  »  f  n 

-  physical  components  of  body  force  vector  (units  of 
force  per  unit  area) 

N^,ig/,N^^ 

-  force-resultants  referred  to  cartesian  parameters 
(units  of  force  per  unit  length) 

-  moment  resultants  referred  to  cartesian  parajneters 
(units  of  force-length  per  unit  length) 

yQ^ 

-  transverse  shear  resultant  referred  to  cartesian 
parameters  (units  of  force  per  unit  length) 

Ti  > 

-  bounding  latitudes  of  paraboloidal  shell 

(") 

-  signifies  quantity  under  bar  is  a  vector 

•d  'm 

-  vertical  slash  before  subscript  signifies  covariant 
differentiation  with  respect  to  metric  of  deformed  body 

xii 


or 


^(X 


n 


m 

f 


1 


comma  "before  su'bscript  signifies  covariant  differen¬ 
tiation  with  respect  to  metric  of  undeformed  body 

Christoffel  symbols  of  second  kind  of  undeformed  "body 


Christoffel  symbols  of  second  kind  of  the  middle  surface 
of  undeformed  body 

Christoffel  symbols  of  second  kind  of  deformed  "body 

Christoffel  symbols  of  second  kind  of  the  middle 
surface  of  deformed  body 

dot  between  two  vectors  signifies  scalar  product 


X 


(ir)J 

n 


cross  between  two  vectors  signifies  vector  product 

second  fundamental  tensors  of  deformed  middle  surface 

second  fundamental  tensors  of  undeformed  middle  surface 
th 

n  Fournier  components  of  displacement  of  V.^,  and 
w^  respectively 

Fourier  integer  in  sin  n  0  ,cos  n  0 

Differential  operators  associated  with  momentless 
behavior  (  i,J  *  1^2^ 5) 


Differential  operators  associated  with  edge  zone 
behavior  (i,j  =  1,  2,  3) 


^  )  -  n^^  Fourier  expansion  for  surface  loads-^^  {  ^ 


a.  , 

1.1 


a 


0 

ij 


k 


4 


12 


coefficients  of  differential  operators 


extensional  stiffness 
bending  stiffness 


X  (R)^ 


large  parameter 


Auxiliary  functions 


xiii 


4 

Roots  of  characteristic  equation;  i  =  1,  2,  --8 

"i  ■ 

Fourth  roots  of  -Ij  i  =  1,  2,  3^  ^ 

m  — 

- 7- 

y  2  -TJ^  ~  inconvenient  constant 

r  = 

km  rjT  - > 

/  yl  +  Y^  -  exponent  of  e 

c . 

1 

Constants  (real)  of  integration  i  =  i,  2,  3^ —  ^ 

Ei(ir)  - 

Auxiliary  functions  (i  =  1,  ..,10) 

i-idf)  - 

Auxiliary  functions  (i  =  1,  2,  3) 

Superscript  e  signifies  edge  zone  portion  of  solution 
Subscript  0  signifies  axi- symmetric  solution  (n  =  q). 
Subscript  1  signifies  anti -symmetric  solution  (n  ■  l) 
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8.1  INTRODUCTION 


Since  an  exact  solution  to  the  set  of  shell  equations 
developed  in  Chapters  I  through  V  [^2]  does  not  seem  possible^  re¬ 
course  must  be  taken  to  solutions  which  will  adequately  describe  the 
behavior  of  a  shell  under  suitably  restricted  conditions.  In 
Chapter  VI  [55], the  shell  was  assumed  to  carry  the  applied  loads 
by  a  momentless  or  membrane  type  of  action.  In  Chapter  VII  [5^]> 
the  shallow  shell  theory,  while  including  the  bending  actions^  restricts 
application  to  shells  whose  middle  surface  slope  is  small.  Unlike 
these  earlier  efforts^  the  present  chapter  will  attack  the  entire 
set  of  shell  equations.  Thus^  although  the  solution  to  be  obtained 
is  approximate;  the  restrictions  inherent  in  the  membrane  and  shallow 
shell  analyses  have  been  removed. 

Let  us  write  the  equations  of  static  equilibrium  of 
forces  (4.4.1.14);  (4.4.1.15)  and  (4.4.1.16)  in  terms  of  the  middle 
surface  displacements.  For  a  complete  sheU  of  revolution;  all 
quantities  involved  must  be  periodic  in  the  circumferential  direction; 
and  Fourier  series  expansions  of  these  quantities  in  0  are  in  order. 
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(sin  ne'j 

|cos  nej 
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w  = 


z 


OO 


z: 


rC05  ne ) 


8.1.1 

8.1.2 

8.1.3 


The  set  of  three  partial  differential  equations  of  force  equilibrium 
can  then  he  reduced  to  a  system  of  Jtn  ordinary  differential  equations. 
For  the  n^^  harmonic  of  the  expansions,  the  equilibrium  equations 


1 


read  as  follows. 


12  ^  ©n^ 


15 


.i-lL:,(v>L;(a,>L:,(w„)l}=^)p®'"^ 
[l,^  (^>l:,(s„)-l;(w„)|-  ^’p,„ 

2^(Upr))+ 


where 


^[LJ;(u.JtLj'^(iiJtL^3(w^)^ 

,4-  /  Z:f  \  '^ 

i  M2(f) 

Pi„(»' -4f'8yr+'?^  ^;„w 


( 1  - 1)^)  p  3.1.6 

Eh  ■‘•sn 


8.1.7 


8.1.8 


^  coefficient  of  the  Fourier  expansion  for 
the  surface  loads  i.e., 


(Sin  ne 
,cos  no 

C05  n© 
5 in  n© 

Sin  ne 
cos  ne 


The  linear  operators  L.  and  L. .  are 
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ij 


1  ~  c 

^11  Si  dif 


'n  dy 


10 


■12 


°  d 
^1^77 


bh-nr  +  b 


0 

10 

0 


I  ^  0  d  ,  ^ 

^15 "  ^11  dy  ‘^fo 

I  ^  _d^  d 

*-  11  dlf^  ^11  ^  ^ 

"  ^11  i"*"  >^10 

1  ->«■  _  d^  d^ 

*-13  ^13773+  ^12  dP 


10 


^11  d^^  ‘^lo 
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o 


+  r 

dil  2£) 


L,,=  t 


“  I  0  d 

J  1-  b 


^zz 


21  dj  ^20 


I  = 

^  25  ^20 


d 


^  21  ^21  ^;^0 
I  ^  =r  L  d  ^  I  d  I 

'-  zz  '^ZZJP  *  °zfgT  *  °zo 

d  ^  d 

^  23  ^2.2  d!^  ^  *^21  dy  '*'  ^20 

L  =c  '  ^  +  c  ° 

^  oi  ^  51  da  ^0 


1  =  b  ° 

^  3Z  ^30 


and 


I  = 

^53  ”  30 

L  *  =  C  —  +  C 

>-  31  '-33  dj3 


d 


32  dj*  ■''^3'd*  ^  *^30 


L,*  '  b 


-i!  +1.  — 

32  d#*-  31  df 


^  'b 


30 


L 


35 


j4-  13  j 

■*■  ^T-2~r~zz  +  0- 


■35  dr 


■52 


dy 


2  1 

+  d^^-rr  +  Q 


■31  dy 


30 


The  coefficients  a..,  w  c..,  a?.^  c?,  (i.  j  »  1, 

ly  ±y  ly  xy  ij'  ij  '  ^ 


2  and  3)  of  the  above  linear  differential  operators  are  all  analytic 
functions  of  and  contain  in  them  the  parameter  n  from  the 

Fourier  expansions.  In  what  follows,  attention  will  he  restricted 
to  n«k.  Only  those  coefficients  needed  in  the  su'bsequent  develop¬ 
ment  are  given  here.  The  remaining  ones  can  he  found  in  the  appendix. 
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-y  + 

(1  +  y^)  V2- 

+  71  (i+v-t-  yd 

~!^(:i(i+v)(H^y+-i^) 

;i(l~2^V 

a7)f^)^/^ 

-y 

(1+y-^) 


>f(]-v) 

2Vf^- 


1 

y__ 

Vi+y^ 

±  n  (5-v) 

2y 

n  (hu) 

2 

(l7V+}f^) 

y(i  +  i^+'%^y  ) 

(TTyT^ 


It  is  seen  that  the  left  hand  sides  of  equations  (8.IA) , 


(8.1.5)  and  (8.1.6)  are  divided  into  two  groups  of  terms.  One  group 
which  is  multiplied  by  k  is  associated  with  the  bending  behavior 
while  the  other  is  associated  with  momentless  or  membrane  behavior. 
The  fourth  power  of  the  parameter  k  is  in  reality  proportional  to 
the  ratio  of  extensional  stiffness  to  bending  stiffness  multiplied 
by  a  factor  porportional  to  the  square  of  a  representative  magnitude ^ 


R,  of  the  radii  of  curvature^  i.e.^ 
If 


h  (X 


extensional  stiffness 


(  R  )‘ 


bending  stiffness 
In  a  homogeneous  shell  of  thickness  h,  the  ratio  of  the  stiffnesses 
is  inversely  proportional  to  the  square  of  the  thickness  while  the 
radii  of  curvature  for  a  paraboloid  is  proportional  to  2f.  For  a 
thin  shell  the  parameter  k  is  therefore  larger  than  unity  since  by 
definition  (see  Section  4.2)  the  ratio  of  2f/h  is  larger  than  10. 

In  general,  large  antennae  have  ratios  of  2f/h  larger  than  200  which 
leads  to  a  value  of  k  greater  than  25. 

The  magnitude  of  k  (for  sufficiently  thin  shells)  and  the 
form  of  the  differential  equations  makes  appropriate  an  asymptotic 
analysis  of  the  solution  to  the  differential  equations.  The  asymptotic 
behavior  of  a  linear  system  such  as  given  by  equations  (8.I.I), 

(8.1.2),  and  (8.1. 5)  has  been  investigated  by  H.  L.  Turrittin  [35] 
and  the  results  of  his  study  form  the  basis  of  the  present  investigation. 
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8.2  THE  AS'jMPTOTIC  BEHAVIOR  OF  A  LINEAR  SYSTEM  CONI’AINING  A 
LARGE  PARAMETER 


Consider  a  real  function  g(k)  of  a  real  variable  k  defined 

over  some  interval  Includlrig  infinity.  A  fonaal  power  series  in  k  ^ 

(denoted  by  P  (k)  ) 

6 


L^O 


with  partial  sum  (k) 

fn 


8.2.1 


•  •)  8.2.2 


is  said  to  be  an  asymptotic  expansion  for  g  (in  the  given  Interval 

of  k)  as  1  k  1—^  oC  (denoted  by  g  "^P^)  if  for  every  m  =  0,  1,  2,  - 

there  exists  some  constant  K  such  tliat 

m 


1 ^ 


J<irL 


8.2.3 


as  |kl->oo  uniformly.  Note  that  if  the  series  P  converges  to  g, 
then  it  is  the  us'oal  Laurent  series  representation  of  g;  but  in 
general  P  (k)  is  not  necessarily  a  convergent  series- 

s 


If 


and  1  f 


6(k) 


g^(k) 

gpCk) 


oC 

g(k)  ~  P^(k)  =  2 


8.2,1* 

8.2.3 


then  we  write 

B.2.6 


whether  or  not  g^  and  g^  have  their  own  asymptotic  expansions.  In 
the  event  that  (and  therefore  g^)  does  not  Ijave  an  asymptotic 
expansion,  (8.2.6)  is  to  be  interpreted  as 


6 


g(k) 


^  P  (k) 
g 


It  is  not  difficult  to  see  from  the  preceding  definition  that 
if  g  has  an  asymptotic  expansion,  then  the  expansion  is  unique.  For 
if  g  has  two  different  asymptotic  expansions  P^^Ck)  and  PgCk) , 


PH)"  Z  -f  .H, 

1  l-mQ  I  H 

CO 


-t 


-L 


8.2.7 

8.2.8 


then 


■*■■■■  8-2-9 


where  and  are  the  partial  sums  of  P^^  and  respectively. 
But  hy  (8.2.5) ^  ve  have 

H!ns,jh-5,„whiir[is,„-^i.i5,„-^i]‘  ^  s.a.io 

Therewith 

^  1  m  K 


Taking  m  =  0,  1,  2 


(in  that  order)  and  letting  k  go  to  infinity 


for  each  m,  we  get  p^  ■  p^^^^  for  all  m. 

On  the  other  hand,  two  functions  g^^  and  gg  which  are  not 
identical  in  a  given  interval  of  k  may  have  the  same  asymptotic 
expansion.  For  instance,  if  equation  (8.2.5)  holds,  then 


“  k 

g(k)  +  e  P  (k)  since 

s 


I,  1  fl  1 ^ 

J^\rr\  l"^!  6  0 

for  m  =  0,  1,  2,  5,  - . 


8.2.15 
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It  can  also  "be  shown  that  (cf.  [36],  ) 

(1)  If  (k) ,  i  =  1,  2,  ni  are  functions  of  k  which 

all  have  asymptotic  expansions,  "the  same 

interval  of  k,  and  f^  - ,  x^)  is  a  polynomial 

in  m  variables,  then  F(k),  defined  hy 

F(k)  -  fp  (g^  (k), - ,  g^(k)  )  8.2.14 

i.e.,  replacing  x^^^  in  f^  by  g^(k)  has  an  asymptotic 
expansion  in  the  given  interval  and  is  calculated  as 
if  all  the  expansions  were  convergent  series. 


(2)  If  g  (k)  is  continuous  for  k  ^  k^  and  possesses  an 


asymptotic  expansion, 

00 


8.2.15 


then 


,  S.2. 
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L-t  i 

(5)  If  the  derivative  of  g  exists  and  is  continuous  and 
if  g'  has  an  asymptotic  expansion,  equation  8.2.15, 
then 


OO 


8-2.17 

In  the  subsequent  analysis^  these  results  will  be  used  frequently  and 
tacitly. 

Next^  let  us  consider  the  second  order  differential  equation 


8.2.18 


where  C|)  (t)  is  either  positive  or  negative  for  all  t,  a  t  —  b. 


and 


8 


k  is  a  large  parameter  such  that 

^(t)| 


8.2.19 


for  all  t  in  the  given  interval.  Intuitively,  the  first  two  terms  in 
(8.2.18)  constitute  the  most  significant  portion  of  the  differential 
equation  as  the  third  term  is  negligible  by  definition.  The  above 
observation  naturally  suggests  a  solution  of  the  form 


where 


x(t)^e' . . Z  Y-i 


f(t)- 


8.2.20 


8.2.21 


and  Y,(t)  can  be  computed  from  a  recurrence  relation  derivable  from 
J 

(8.2*l8)  and  (8*2,20).  Since  the  series  on  the  right  hand  side  of 
(8.2.20)  is  not  necessarily  a  convergent  series  in  k  for  a  ^  t 
it  is  therefore  (formally)  an  asymptotic  representation  of  x(t).  For 
discussions  of  the  asymptotic  behavior  of  (8.2.18) ^  readers  are  referred 
to  [58],  [39].  [40]  and  [4l]. 

As  a  natxiral  extension  of  the  results  obtained  for  the  single 
ordinary  differential  equation  (8.2.18) ^  Turrittin  investigated  the 
asymptotic  behavior  of  a  set  of  simultaneous  linear  (real)  ordinary 
differential  equations 

X  =  ^  °  A(t,'^)x  8.2.22 

ot 

where  m^5;  1  is  an  integer;  x'  and  x  are  m  x  1  column  matrices;  A  is 
an  m  X  m  square  matrix;  k  is  a  real  nmber  which  is  large  relative  to 
imity.  Each  element  of  the  matrix  A  (t,k)  is  assumed  to 

have  an  asymptotic  expansion  in  k 
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8.2.25 


oo 


V'  I  -i  i 

.(t)^  L  -h  (X  .  .  (i) 

i-.o  ij 


I 


where  the  cK.  j  (■t)'s  are  infinitely  differentiable  in  t.  Denote  the 
matrix  whose  elements  are  (X  >  i*e., 


0  L  . 


8.2.24 


and  the  characteristic  roots  of  A  (t)  hy  An  (t),  q  =  1,  2,  m. 

o  f 

For  a  system  for  which  A^(t)  has  the  same  number  of  distinct  charac¬ 
teristic  roots  for  all  t,  a  4  t  ^b  (i.e.,  if  A^(t)  and  Xj(t)  are  two 


roots,  then  either 


or 


A,(t)=  A  At) 
^  i 


A^(i)  ^  A^  (t) 


8.2.25 


8.2.26 


for  all  t,  1)  4  t  a) .  Turrittin  showed  that  there  exists  a  solution 
is  a  column  matrix)  of  the  system  (3,2.22)  such  that 
each  of  its  elements^  asymptotic  expansion  of  the  form 


I  p,o  '  tr 


_|n 


8.2.27 


1^- 


8.2.28 


for  jk|  ^1,  where 

^ir^  (t)  and  q^'^^(t)  are  computed  by  a  certain  algorithm.  As  j  varies 
from  1  to  m  (to  account  for  all  the  m  characteristic  roots),  we  will 
have  the  desired  asymptotic  expansions  of  the  m  linearly  independent 
solutions  to  the  given  system  of  differential  equations. 
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8.5  the  eigenvalues  of  A^(t) 

As  the  first  step  in  the  application  of  the  results  of 
Section  8.2  to  our  shell  problem^  we  must  find  the  characteristic  roots 
of  the  matrix  corresponding  to  A^  forthe  system  of  equations  8.1.4,  8.1. 5> 
and  8.1.6.  In  order  to  do  this,  we  will  convert  our  system  of  differ¬ 
ential  equations  into  the  form  of  (8.2.18).  With  no  loss  in  generality, 
we  may  write 


8.5.1 

11 

C 

CD 

(k 

8.5.2 

Wn 

8.5.5 

dlf 

8.5.4 

d  K 

11 

8.5.5 

dw,-, 

d!( 

8.5.6 

8.5.7 

- 

8.5.8 

where  the  k  in  these  equations  is  the  shell  parameter  defined  by 
equation  (8.1.7)^  and  the  take  on  the  role  played  by  in 

section  8.2.  We  then  have  the  following  system  of  homogeneous  first 
order  differential  equations  which  is  equivalent  to  the  set  of  shell 
equations. 

®-5-5 


11 


dr 

- 

"dv 

A§^= 

dlf 

dir 

= 

dr 

dr 

dr 


more  concisely^ 

M. 

dr 

vhere 

0 


A^(t,k)-  0 


0 

0 

0 


5 


A  fS, 


8.5.10 

8.5.11 

8.5.12 

8.5.15 

8.5.IU 

8.5.15 

8.5.16 


8.5.17 


0 

0 

0 

0 


0 

0 

0 


0 

0 

0 

0 


20 


22 


0 

0 

-a! 


JO 


'54 


1 

0 

0 

0 

-c' 


21 


22 


0 

0 


^5U 


0 

1 

0 

0 


0 

0 


-a 


0 

0 

1 

o 

11 


"12 

0 

0 

0 


0 

0 

0 

0 


1 

0 


8.5.18 


8.5.19 
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The  eight  characteristic  roots  of  this  matrix  are 


X^(^) -  \^{x) -  0  8.5*20 

and  K(,Cf),  A^(v),  and  Xg(  S'  )  which  are  the  four  roots  of 

K ‘^51  -  Cl  =  0 

(CIj4Cj2^  d^j  8.5*21 

Note  that  these  roots  satisfy  the  conditions  (8.2.25  and  8.2-26). 
Moreover,  they  are  independent  of  n,  the  parameter  from  the  Fourier 
expansions.  With  these  characteristic  roots,  we  may  proceed  to  compute 
the  asymptotic  solutions  to  our  system  of  differential  equations  in 
accordance  with  an  algorithm  outlined  in  [55]*  Inasmuch  as  this 
algorithm  is  rather  involved,  we  shall  adopt  a  simpler  procedure 
which  is  possible  because  of  the  nature  of  our  system  of  equations. 

Observe  that  for  our  problem,  m^  =  1  (cf.  (8.5«17)  )i  the 
equation  corresponding  to  (8.2.28)  is  therefore  reduced  to 


for  all  j.  With  no  loss  in  generality,  we  set  q^*^^ 
it  does  not  involve  k;  therewith 


8.5*22 

( ^  )  ■  0  because 


da^r  j  . 


8.5.25 

(8.2.27),  and  (8.2.28)  for  the  case  at  hand  then  combine  to  yield 


co 


'  L  '  '  If 


-r 


8.5.24 


/,(/)  (p 

for  i,  J  =  1,  2,  — ,  8  and  take  on  the  role  of  .  In  other 

words,  the  work  of  Turrittin  assures  us  that  there  exist  solutions 
to  (8.1.4),  (8.1.5),  and  (8.1.6)  having  (8.5*24)  as  their  asymptotic 
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expansions.  We  will  see  later  that  these  asymptotic  solutions  divide 
themselves  into  two  groups.  The  leading  terms  of  the  expansions 
associated  with  the  four  zero  roots  (8.5*20)  will  identify  themseives 
with  the  membrane  and  inextensional  solutions  (i.e.,  displacements 
associated  with  the  homogeneous  solutions  to  the  membrane  force 
equilibrium  equations  and  to  the  strain-displacement  equations) .  The 
second  group  of  the  asymptotic  solutions  associated  with  the  non-zero 
roots  (8.5*21)  will  be  shown  to  be  significant  only  in  a  region  near 
the  edge(s)  of  the  shell  since  their  magnitudes  decay  rapidly  with 
distance  from  the  edge(s)  and  is  thus  referred  to  as  the  edge  zone 
solutions.  Consequently,  the  membrane  and  inextensional  solutions 
dominate  the  behavior  of  the  shell  interior  and  together  they  are 
referred  to  as  the  interior  solutions.  We  will  discuss  separately  these 
two  groups  of  solutions  in  the  following  two  sections. 


Ih 


8.4  THE  INTERIOR  SOLUTIONS 

For  A y ) ;  i  =  1;  2,  3  and  4;  i.e.^  the  zero  eigenvalues^ 
(8. 3 -24)  becomes 


^  8.4.1 

f  5  =  0  t 

a.h.S 


and 


8.4.4 


To  find  the  asymptotic  homogeneous  solution  to  the  given  system  of  shell 
equations,  we  substitute  (8.4.2),  (8.4.3)  and  (8.4.4)  into  (8.1.4), 
(8.1.5)  and  (8.1.6)  with  the  forcing  terms  (i.e.,  the  right  hand  side) 
suppressed  and  collect  coefficients  of  the  same  power  of  k.  Equating 
these  coefficients  to  zero,  we  get 


L^,(/S3,)  =  o  8.4.5 

^  m?  ’’’ ^rni  ^^10^  ~ 

^  ni3  mi  ^  (^)l^  ~  ^  8.4.7 

and  for  S  ^  3 

-  {L*3  %-->>  L;  L  *  (,6,5-5)]  8.1..8 

where  m  =  1,  2  and  3*  Quantities  with  a  negative  subscript  should  be 
taken  as  zeros.  In  general,  can  not  satisfy  the  three  equations 
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given  by  (8.4.5)*  Therefore,  we  must  have 


8.4.9 

Similarly,  ^  /^lO  general  can  not  satisfy  the  three  equations 

given  by  (8.4.6);  therefore 

*  ^,o(ir)=  0  8.4.10 

In  view  of  (8.4.8),  (8.4.9)  end  (8.4.10),  we  may  now  (with  no  loss  in 
generality)  rewrite  (8.4.2),  (8.4.5)  end  (8.4.4)  in  the  form  of 


oO 


-5=0 

8.4.11 

OO 

-5*0 

8.4.12 

and 

oO 

s*o 

8.4.15 

so  that 

(8.4.7)  and 

(8.1j 

l■.8)  become 

=  0 

8.4.14 

and 

5 

3 

)  j  0 

8.4.15 

for  m  = 

1,  2  and  5. 

The 

j  set  of  equations 

given  hy  (8. 

4.l4)  upon 

comparison  with  the  complete  set  of  shell  equations  (8.1.4),  (8.1.5) 
and  (8.1.6)  can  be  recognized  as  the  membrane  equations  (written  in 
terms  of  the  middle  surface  displacement  components  by  way  of  the 
strain-displacement  relations)  since  the  terms  containing  k  are  not 
present.  The  four  linearly  independent  solutions  to  these  equations 
(corresponding  to  the  fact  that  zero  is  a  quadruple  root)  are  the  same 
as  those  previously  obtained  in  Chapter  VI.  Since  two  of  these 
linearly  independent  solutions  are  the  homogeneous  solutions  of  the 


l6 


strain-displacement  equations  while  the  other  two  are  associated  with 
the  homogeneous  solutions  to  the  membrane  force  equilibrium  equations ^ 
they  are  often  referred  to  as  the  inextensional  and  membrane  solutions^ 
respectively.  While  in  principle  the  higher  order  terms  in  the  ex¬ 
pansions  can  be  obtained  by  way  of  (8.4.15) ^  the  calculations  are^ in 

/  4 

general^ tedious .  Since  the  additional  terms  are  of  order  1/k  and 
higher ;  we  will  only  be  interested  in  the  leading  term  of  each  expansion 
in  the  subsequent  development^  i.e.,  we  will  only  be  concerned  with 
the  above  described  solutions  to  (8.4.14). 
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8.5  THE  EDGE  ZONE  SOLUTIONS 


For  X  ^ ( 0  ^  i  “  5;  6,  7  and  8,  i.e. , 

the  non- zero  eigenvalues 

(8.5.25),  we  have 

""  r° '  ’f 

8.5.1 

8.5.2 

Ci)  r  (i)  --i 

^  1=0  1  ^ 

8.5.5 

Again  we  substitute  these  expressions  into  (8.1.4),  (8.1. 5)  and  (8.1.6) 
and  collect  the  coefficients  of  the  same  power  of  k.  Equating  each  of 
these  coefficients  to  zero;  we  have  the  following  recurrence  relations: 


8.5.4 

8.5.5 

8.5.6 

8.5.7 

8.5.8 

8.5.9 


18 


19 


1  ^ 


a) 


2.A 


K 


^  K  ^ 

8.5.10 


/ 


w' 


i  '  i-'  ■  ■“'=  -■  I’^'i 

[^13^  jS^. 

8.5.11 


a';/// 


(i)// 


■*■  “•  II  /-3 


A)/  ^  (0 

*  N 


11 


0 


11 

|-4 

10 


b.i  ^ 


(0// 


w/ 


5^. 

1-’ 


li  /^i|  +  <^11  (^1/  ■"  ^'»  f^ij 

(i) 


a) 


A  -2^ 

0  ■  =  a 


0  ^  (0 


A 

^'-'  ■ 

z^> 

V 


^ ^ 

0  ^^  /}(.€)! 


8.5.12 

8.5.15 

8.5.1^ 

8.5.15 


a) 


°  /l  (t)//  I  ° 


j  ^0  f'''-^ 

8.5.16 
6’) 


a'  ^ 

^ZZ  ^  + 


‘^;i2  ‘^Zl  ^ 


0) 


8.5.17 

8.5.18 

I/)  0) 


Z1 


^■kz  *  [^22  t>21  ^ 


8.5.19 
(0 


(0  .  (0/  . 
'21  f^lj  ^ZO  (^tj 


8.5.20 


20 


8.5.28 


50 

.i-C> 


vhere 


(n 


^0  ;^Aj 


'f  -  f  X  •  dv 


8.5.29 

Quantities  with  negative  subscripts  are  again  to  be  taken  as  zeros. 

For  j  =  0,  equations  (8.5.7)>  (8'5.8)^  and  (8.5.9)  become 


8.5.30 


and 


bj)  (S^c  *  ^2°  (^.»)  ■  ° 

8.5.32 


1^.0 

respectively.  (8.5.30)^  (8.5.31) ^  and  (8.5.32)  have  nontrivial  solutions 
for  jS  2_o>  ^20'  (^30  determinant  of  the  coefficients 

vanishes : 

a.Ya'(^)- 

t  \  L 

vhich  is  assured  by  (8.3*21) 


o  ^0 
^31  ^11 


a 


a  r  ^ 

^34  ^  12- 


a 


34  J 


=  0 


8.5.33 


For  j  “  1,  (8.5.7) ^  (8.5.8),  and  (8.5.9)  become 


“■1°  A,  <".1  A, 


50 


^  0 
+  a 


10  50 


2X-  C. 


0  \  ®  \  ^ 

c  A.  +  c  A 

1 10 

L  11  L  >2  t  J 

0 


8.5.34 


21 


1^31  ^  ^L  f>l 


andp 

t 


/^3.^C5>c^n"  ^°-^'^\§ 


/ 

30 


{4a „  A.'  A^ a„  A=  ]  +  c;  ^,'  ^  C,:  /5,„  =  0  8.5.36 


where 


( 

dif 


The  first  and  third  equations  yield  (keeping  in  mind  the  results  for 


j  =  0) 


/«»(») 


=  exp^'— / 


+  ^35  -  (‘^.^il  *^10  ^11  *“30  ) 


X; 


a 


44 


Cu  x^ 


.  c,-  a,r 

('^il)^  ^34 


S;^34V 


dY 


8.5.37 


From  {8.5.30)  and  (8.5.31)^  we  have 


/S,o(«)  = 


=  - 


30 


8.5.58 


and 


b;,  X;"  (C  A»)  ®-5-39 

Terms  associated  with  higher  powers  of  (k)~^  to  each  expansion  can 
again  he  computed  hy  way  of  the  recurrence  relations  (8.5.^)  through 
{8.5.6)^  although  the  computations  are  not  only  tedious  hut  may  he 
very  difficult  because  new  sets  of  differential  equations  arise.  Since 
k  is  large,  these  higher  order  terms  are  negligibly  small  in  comparison 
with  the  leading  term  and  therefore  will  not  he  considered  in  the 
subsequent  development,  i.e,,  we  will  consider  only  ^20’  ^^^^50 
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8.6  THE  PARTICULAR  SOLUTION 


For  a  shell  under  gravity  loading,  the  components  of  the  load 
intensity  vector  are  [cf.  [55I  ) 


K 


Vi+r 


sin  Y  sine  -!(  C05'^ 


8.6.1 


and 


^  If  sin  sm  e  C05  Y 


8,6.2 


8.6.5 


The  form  of  the  corresponding  )  (  i  ■  2,  5)  suggests  that 

any  particvilar  integral  of  the  dependent  variables  cannot  contain  a 

^  ) 

factor  e  '  for  any  X  ( )  which  is  not  identically  zero.  Thus  we 

will  seek  particxilar  asymptotic  solutions  of  the  form 


rnp  '  %  »'} 

(n  =  0,  1) 

8.6.4 

00 

(n  =  1) 

8.6.5 

(n  =  0,  1) 

8.6.6 

Substituting  (8.6.4),  (8.6.5),  and  (8.6.6) 

into  (8.1.4), 

(8.1.5),  and 

(8.1.6),  collecting  coefficients  of  the  same  power  of  k  and  equating 
them  to  zero,  we  get  for  n  -  0 


4f 

Eh 

“Eh 


cos  Y 


cos  Y 


If 


8.6.7 


8.6.8 


25 


and 


(53=  I)  8.6.9 


where  m  =  1  and  Q-nd  for  n  =  1 

I  f'  \  \  f-  \  \  f-  '  -  4f^(l-z;*)Al^5'nV  r 


■|5  150'  12  120'  ‘“11'  iio' 


 4f^(l-y"V^5in^ 


8.6.10 


Eh 


'/O-hSinY  r  - 1 

^ - [yyvK^J 


L  (S  VL  «  )+L  )-  ^  M 

35  150  52  120'  3r  HO^  ^ 


th 


8.6.11 


8.6.12 


and 

8-<5-13 

where  m  =  1^  2^  and  5* 

Clearly^  the  leading  term  of  the  expansions  given  by  (8.6,4)^ 

(8.6.5) >  and  (8.6.6)  can  be  identified  as  the  particular  solution  to 

the  corresponding  middle  surface  displacement  in  the  membrane  analysis 

of  a  shell  loaded  by  gravity.  Higher  order  terms  of  these  expansions 

can  be  obtained  from  (8.6.9)  Q^nd  (8.6.13).  However^  since  these 

1 

additional  terms  are  of  order  — or  higher^  ( small  in  comparison 

Ic 

vith  the  leading  term)  they  are  therefore  omitted  from  further  con¬ 


2k 


siderations. 


8.7  SOME  RELEVANT  RELATIONS 


Some  of  the  relations  derived  in  Chapters  I  through  V  will 
he  recapitulated  here  in  a  form  suitable  for  the  purpose  of  asymptotic 
solutions . 


4.  J2L 

dir 


V  V 

}  8.7.1 

ae  J 


N  V  ^  4 

^  (1-D^)  I.  L  2f  (i-(-5^)V^ 


1 

2/yi7s^ 


w  + 


^  .  1  ^6  ) 
di'  ^  de  J 


zjnYWP- 

8.7.2 


=  E  In  C  1  dltp  _ _ Hft  ^ _ 1  due } 

•"6  2(1 +d)  izfi  de  zjii^  J 

8  •  Y  •  5 


=  -b  [  !  ^ 

4f^l('i+if^)  d(f^ 

4.  IL.  h^.  4  1  4.  ^  ] 

yYhT^  So  (1+y^)^  dJ  y(i+y^)^ 

8.7.4 

=  -b  [  XI  a^w  f  d^w  [l  +  ^j-x) dv/ 

&  dV'^  f  dd""  yO+y^)^  dY 

^ _ l_  ^  V  dur  ^  ^  1 

yVITy^  de  dif  Y(ihy^)^  "'J 


8.7.5 

=  -b(i-x)  (  1  dV  _  1  dw  1  dcte 

4p  dYde  H'Vl+Y^  de  “2XFFJ  dY 

_  ,  1  du^) 

2Y(ltY^)^  ®  2!f(l+Y^)^^  de  J 


8.7.6 


25 


Q  =  1  (  ^(2;fyMre)  ^  ^  of  ^ 


re 


8.7.7 


n  1  (  ()(^f^^Mr)  .  Vi^Mre)  .01^} 

4f^jyFFl  ^ 


LL, 


1  (3w 


_  r  u-ft  + _ L.  ^ 

^e~  |2fYTt7^  2/K  ae. 


8.7.8 


8.7.9 


8.7.10 


where 


D  = 


Eh- 


U(l-TJ^) 


8.7.11 
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8.8  THE  LEADING  TERM  OF  ’THE  ASYMPTOTIC  EXPANSION  FOR  THE  EDGE  ZONE 
SOLUTIONS 


We  shall  consider  in  this  section  asymptotic  solutions 
for  a  paraboloidal  shell  of  revolution  corresponding  to  the  four  charac¬ 
teristic  roots  X  •  (  ^^  )  which  are  not  identically  zero. 


where 


X- (x)  dx  =  xj.  m  f 


dx 


8.8.1 


8.8.2 


and  (J-  is  the  appropriate  root  of  the  equation 
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For  the  n  harmonic  of  the  Fourier  expansion^ 
the  expansions  (8.5.I);  (8.5*2)  and  8.5.5) 


the  leading  terms  of 
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where  C  is  a  complex  constant  of  integration.  Retaining  only  the  leading 
terms  of  the  asymptotic  solutions^  the  relevant  edge  zone  stress  resultants 
and  bending  moments  can  be  obtained  from  the  middle  surface  displacements 
by  way  of  the  relations  listed  in  section  8.6.  In  what  follows^  the 
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superscript  e  attached  to  a  dependent  variable  means  the  edge  zone  part 
of  that  dependent  variable.  Similarly  the  superscript  i  refers  to  the 


interior. 
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i  =  5^  C,  7;  8  constants  of  integrations. 

e  “1 

The  leading  terra  of  is  proportional  to  k  .  Since  we 

0 

shall  have  no  need  for  it^  the  explicit  expression  for  will  not  he 
given  here.  To  compute  the  coefficient  associated  with  this  leading 
terra  would  necessitate  the  derivation  of  the  second  term  of  the  ex¬ 
pansions  (8.5.1)  and  (8.5.3) •  On  the  other  hand;  the  error  accrued  to 
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the  neglect  of  is  of  the  order  (k"^)  in  the  presence  of  the 
interior  solution  corresponding  to  =  0.  Observe  that  the  effects 
of  the  solutions  obtained  in  this  section  due  to  the  presence  of 
<S  are  of  significance  only  in  a  region  near  the  edges  of  the 

shell;  hence ^  they  are  often  referred  to  as  the  edge  zone  solutions. 
In  many  situations^  the  edge  zone  is  only  a  small  percentage  of  the 
span  of  the  shell. 
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8.9  THE  LEADING  TERMS  OF  THE  ASYMPTOTIC  SOLUTIONS  FOR  AXI-SYMMETRIC 
GRAVITY  LOADING 

In  this  section  the  interior  solution  (section  8.4  and 
chapter  VI),  the  edge  zone  solution  (section  8. 8), and  the  particular 
solution  (section  8.6  and  chapter  Vl)  will  be  combined  for  a  paraboloidal 
shell  under  axi- symmetric  gravity  loading.  It  should  be  recognized  that 
the  result  represents  only  the  leading  terms  (section  8.4  and  8.5)  of 


the  asymptotic  solution  to  the  shell  equations  8.1.4,  8.1.5^  and  8.1.6. 
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In  addition,  ve  will  also  be  interested  in  the  rotation  60 
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and  have  been  omitted  in  the  expressions  (8.9.2) 

and  (8^. 9.5)  for  and  N^,  respectively,  because  they  are  of  the  order 

k  ^  in  the  presence  of  the  interior  solution  (see  equations  (8.8.8)  and 

(8.8.15)  ).  In  contrast,  60^,  ,  and  which  are  of  order  k 

r  r  Q  r 

-2  -2  -5 

k  ,  k  ^  and  k  respectively,  relative  to  the  edge  zone  solutions, 
cannot  be  neglected  on  the  relative  order  of  magnitude  basis  because 
they  are  the  dominant  terms  in  the  interior  of  the  shell.  Nevertheless, 
these  same  terras  which  represent  the  bending  actions  of  the  shell  have 
been  omitted  since,  relative  to  the  membrane  actions,  they  do  not  con¬ 
tribute  significantly  to  the  stresses  and  distortions  of  the  shell.  An 
alternate  derivation  of  the  above  results  can  be  found  in  [12]. 


The  six  real  constants  of  integration  C^,  Cg,  C^,  Cg, 
and  Cg  can  be  chosen  in  such  a  way  that  an  admissible  set  of  edge 
conditions  is  satisfied.  Considered  herein  are  boundary  conditions 
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which  are  commonly  called  (a)  simply  supported,  (b)  clamped,  and  (c)  free. 

The  simply- supported  edge  at  a  fixed  value  of  (i  requires 


the  following: 
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The  clamped  edge  at  a  fixed  value  of  X  requires  the 
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The  free  edge  at  a  fixed  value  of  V  requires  the 
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Fccr  a  more  detailed  explanation  of  the  houndary  conditions, 

see  Section 
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8.10  THE  LEADING  TERMS  OF  THE  ASYMPTOTIC  SOLUTIONS  FOR  ASYMMETRIC 
GRAVITY  LOADING  (  T  =  ^  ) 

The  leading  terms  of  the  asymptotic  solutions  of  a  shell 
whose  axis  of  revolution  is  perpendicular  to  the  gravity  vector  arc 
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In  the  derivation  of  the  above  results,  terms  of  the 
order  (k)  ^  m  ^  1  are  again  neglected  in  the  presence  of  unity  for  the 
various  reasons  discussed  in  section  (8.9)*  The  eight  constants  of 
integration  C^,  i  =  1,  -  -  -  8,  are  to  be  determined  by  the  boundary 

conditions  which  are  commonly  called  (a)  simply  supported,  (b)  clamped, 
and  ( c)  free • 

The  simply- supported  edge  at  a  fixed  value  of  ^  requires 
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The  clamped  edge  at  a  fixed  value  of  Y  requires  the 


following: 
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8.11  THE  LEADING  TERMS  OF  THE  ASYMPTOTIC  SOLUTIONS  FOR  A  SHELL  UNDER 
UNIFORM  PRESSURE 


The  membrane  behavior  of  a  paraboloidal  shell  of 
revolution  under  uniform  pressure  has  been  investigated  in  [^3^*  The 
results  presented  there  coupled  with  the  results  given  in  section  8.7 
(with  n  »  0)  constitute  the  leading  terms  of  the  asymptotic  expansions 
of  the  stress  resultants^  bending  moments,  and  raid- surface  displacements 
which  describe  the  shell  behavior  under  uniform  pressure.  If  the 
intensity  of  the  pressure  is  with  respect  to  the  raid- surface  of 

a  shell  whose  normal  is  assigned  to  be  positive  when  directed  inward, 
we  have 
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where  T  ,  k,  and  m  are  defined  hy  (8.7.20) ^  (8.1.25) ,  and  (8.7*2), 
respectively.  C i  =  1,  2,  - 8,  are  constants  of  integration. 
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8.12  THE  SINGUIARITY  AT  THE  APEX 


The  presence  of  a  large  parameter  k  has  led  to  an 

asymptotic  analysis  of  the  solution  to  the  shell  equations.  Implicit 

in  this  analysis  is  an  assumption  concerning  the  order  of  magnitude  of 

the  coefficients  associated  with  the  differential  equations.  Observe 

2 

■that  if,  in  (8.2.18),  <|>  (t)  vanishes  for  some  t^,  a  ^  h,  while 

2  2  2 

(^0(^0)  ^0^  then  whatever  the  magnitude  of  k  ,  the  term  k  <p  fails 

2 

to  remain  significant  relative  to  ^  at  t  •  t^.  In  fact,  if  4^  is 
continuous,  then  for  a  fixed  finite  k,  however  large,  there  exists  a 
S  such  that 


^  I  ^(t) 


8.12.1 

for  |t“tol  S  .  Thus,  for  lt“tol  <5  ,  the  formal  asymptotic  solution 

(8.2.20)  would  cease  to  he  meaningful. 

An  examination  of  (8.1.^)-),  (8.1. 5),  and  (8.1.6)  reveals 

that  some  of  the  ^  ^ij^  ^ij  =  0,  i.e.,  the  apex  of 

the  shell,  while  the  corresponding  a..,  h  ,  and  c.,  do  not.  Thus, 

ij  ij  ij 

there  will  he  some  neighborhood  surrounding  the  apex  in  which  the 
asymptotic  solutions  obtained  in  section  8.8  cannot  be  used  to  describe 
the  behavior  of  the  shell.  A  method  of  asymptotic  integration  has  been 
developed  by  R.  E.  Langer,  yielding  a  solution  which  is  valid  at  the 
singularity  K  =  0  (cf.  [4o],  [4l]  ).  Using  Langer 's  results,  C.  N. 
DeSilva  investigated  the  axi- symmetric  deformation  (including  the  trans¬ 
verse  shear  deformations)  of  a  paraboloidal  shell  of  revolution  and 
derived  the  leading  terms  of  the  asymptotic  representations  of  the 
rotation,  63^(5’  ),  and  the  horizontal  stress  resultant,  H  (  ^^  ) ,  which 
are  valid  for  0  ^  a^y^b-oO[ 42] .  The  other  stress  resultants,  bending 
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mornents  and  middle  surface  displacements  are  derivable  from  and  H. 
These  quantities  are  all  in  the  form  of  the  appropriate  Bessel  functions 
with  argument  k  A  ( if  )  where  0  as  S'  0  .  It  is  not  difficult  to 

verify  that  the  results  given  in  section  8.8  for  n  =  0  are  just  those 
obtained  by  DeSilva  with  the  Bessel  functions  replaced  by  the  first  term 
of  their  asymptotic  expansions  and  the  transverse  shear  deformability 
suppressed.  While  the  asymptotic  representations  of  the  Bessel  functions 
cease  to  be  meaningful  as  Y  ^  0  ^  the  Bessel  functions  themselves  remain 
unaffected. 

An  asymptotic  solution  to  the  shell  equations  which  is 
valid  at  the  apex  and  without  the  restriction  of  axi- symmetry  has  also 
been  obtained  by  C.  R.  Steele  (cf.  [43]  ).  But^  since  the  region  in 
which  our  asymptotic  solutions  (given  in  section  8.8)  are  erroneous  is 
extremely  small  for  a  thin  shell  (see  discussions  of  results)  and  since 
the  exact  behavior  of  the  shell  at  the  apex  can  be  deduced  from  the 
differential  equations^  the  behavior  near  the  apex  can  be  obtained  to 
a  good  approximation  by  simple  interpolations.  For  engineering  purposes^ 
it  is  not  necessary  to  perform  the  proposed  refined  asymptotic  analysis^ 
although  such  an  analysis  would  offer  a  better  understanding  of  the 
behavior  of  the  shell  near  and  atihe  apex.  For  instance^  it  can  be  con¬ 
cluded  from  the  results  of  this  refined  analysis  that  the  condition  of 
finiteness  at  Y  =  0  requires  the  vanishing  of  the  constants  Crj  and  Cg 
in  the  results  listed  in  section  8.7*  The  relatively  large  amount  of 
computer  time  needed  for  the  use  of  the  double  precision  Bessel  function 
subroutines  was  also  a  strong  factor  in  our  decision  not  to  extend  the 
present  analysis. 
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INTRODUCTION  TO  NUMERICAL  RESULTS 


Numerical  calculations  are  presented  in  curve  form  in  the  sub¬ 
sequent  nine  sections.  It  is^  of  course^  impossible  to  anticipate  all  the 
configurations  which  may  be  of  interest  to  the  antenna  designer.  However, 
a  fairly  wide  range  of  sizes  are  covered  and,  additionally,  the  results 
are  presented  in  normalized  form.  These  results  can  be  used  to  determine 
the  magnitude  of  displacements  (and  stresses)  for  purposes  of  preliminary 
design.  A  more  detailed  analysis  of  a  specific  configuration  can  be  made 
by  use  of  the  digital  computer  programs  which  have  been  developed  for  this 
report. 

The  results,  in  addition  to  being  used  for  numerical  data,  also 
illustrate  the  behavior  of  paraboloidal  surface  structures.  Most  of  the 
text  accompanying  each  of  the  nine  groups  discusses  various  aspects  of 
the  shell  behavior  which  have  been  revealed  by  the  numerical  calculations. 
In  general,  the  shells  exhibit  the  same  qualitative  behavior  as  that  given 
by  a  shallow  shell  analysis (reference  Part  III)  but  there  are  some  signifi¬ 
cant  differences. 

An  attempt  has  been  made  to  present  the  results  in  a  uniform 
manner  but  deviations  will  be  noted  to  accommodate  the  specifics  of  each 
group.  Dimensions  and  boundary  conditions  are  presented  in  tables.  Both 
the  symmetric  and  anti-symmetric  behavior  under  gravity  loads  are  shown 
with  the  former  corresponding  to  a  pointing  angle  of  =  0°  and  the 
latter  to  •  90°.  The  curves,  in  addition  to  three  groups  of  numbers, 
bear  the  letter  s,  meaning  symmetric  and  a,  meaning  anti- symmetric,  e.g., 
figure  8.15.1.1s  and  figure  8.15«5-la*  Some  of  the  parameters  such  as 


Q*  and  M*  tend  to  overlap  if  plotted  on  the  same  scale.  These  have  been 
separated  in  the  interest  of  identifying  each  curve  but  there  may  be  some 
confusion  regarding  the  scales.  However,  a  careful  second  examination  of 
each  figure  should  be  sufficient  to  relate  each  curve  with  the  proper 
ordinate  and  abscissa  scales. 


The  normalizing  factors  are  as  follows  (see  section  Y.12): 


8.15  Group  I  --  Four  Closed  Shells  with  Outer  Edges  of 

V  =  .20,  .UO,  .80,  1.20 

The  Group  I  results  presented  in  this  section  illustrate  the 
behavior  of  four  shells,  closed  at  the  apex,  which  differ  only  in  their 
outer  radii.  Dimensions  and  other  pertinent  data  are  presented  in 
Table  8.I5.I.  The  four  outer  radii  are  250eU,  U60.8,  921.6,  and  I382.U 
inches  which,  since  the  focal  length  is  576  inches,  result  in  slopes  at  the 
outer  edges  of  Y  ®  *20,  .Uo,  .80,  and  1.20,  respectively. 

An  examination  of  figure  8.13.1s  reveals  the  very  interesting 
result  that  the  maximum  deflections  are  a  somewhat  unconventional  function 
of  the  outer  radius.  A  table  of  maximum  deflection  (w*  )  for  different 

shells  under  symmetric  loading  vs.  the  position  where  it  takes  place  (Y^) 
is  shown  below: 


Table  8.I3.2 
Symmetric  Loading 


u  2 

max 

.20 

-.U17 

.144 

■  ho 

-  .411 

GO 

.80 

-  .609 

0 

1.20 

-L043 

0 

We  observe  that  case  2s  with  the  outer  edge  at  62  =  -^0  experiences  a 
w^x  'Which  is  actually  less  than  that  experienced  by  the  smaller  shell 
(case  Is).  Furthermore,  case  ^is.  which  is  a  shell  having  a  diameter  three 
times  that  of  case  2s,  experiences  deflections  only  2.5  times  that  of 
case  2s.  Generally  speaking,  with  the  same  relative  increase  in  length, 
the  transverse  deflections  of  beams  and  plates  which  vary  as  the  square 
of  the  span  length  would  yield  deflections  nine  times  greater.  Hence, 
these  results  clearly  illustrate  the  stiffness  advantage  of  a  shell 
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structure  over  that  of  heams  or  plates. 

Note  also  that  the  edge  zone  hehavior  for  the  two  smaller  shells 
causes  a  hump  in  the  deflection  curve  which  places  the  maximum  deflection 
near  the  outer  "boundary .  The  hump  is  more  pronounced  in  case  Is  than  in 
2s.  The  edge  zone  "behavior  is  also  clearly  visi"ble  for  the  two  larger 
shells  although  the  hump  is  overshadowed  "by  the  interior  "behavior. 

A  "better  indication  of  the  edge  zone  "behavior  is  given  "by  the 
*  *  * 

plots  of  Mg,  ,  and  Q^.  We  o"bserve  that  d,  the  width  of  the  edge  zone, 
is  approximately  the  same  for  all  four  sizes  and  agrees  well  with  the 
result,  ^  .006°  (see  equation  7*12. l.l)  given  "by 

8.13.1 

The  maximum  values  of  these  quantities  which  characterize  the  "bending 

-K*  -X-  *)(* 

"behavior  occurring  in  the  edge  zone,  i.e.,  ,  M  ,  and  Q  ,  decrease 

y  r  r 

as  the  outer  radius  increases.  This  means  that  the  "bending  stresses  will 
be  larger  for  a  smaller  shell - 

y 

The  manner  in  which  the  curves  of  N  and  coincide  should 

r  9 

also  be  noted.  This,  again,  illustrates  the  membrane  behavior  of  the 

shell.  As  has  been  previously  described  (see  sections  7-12.1)  the  shell 

prefers  to  carry  the  applied  loads  by  membrane  action;  bending  action 

takes  place  near  the  edges  in  order  to  satisfy  practical  geometrical 

constraints.  Hence,  the  interiors  of  shells  will  develop  nearly  the 
*  * 

same  values  of  N  and  providing  the  edge  zone  widths  are  not  a  major 
r  y 

proportion  of  the  total  radii. 

The  anti- symmetric  results  are  shown  in  figures  8.13.5  through 
8.13.11.  We  note  for  all  except  the  smallest  shell  that  the  displacements 
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for  the  anti -symmetric  gravity  loading  are  larger  than  for  the  symmetric 
situation.  For  the  two  larger  shells ,  the  maximum  anti- symmetric  deflections 
are  more  than  twice  as  large.  This  can  he  attributed  to  the  unfavorable 
orientation  of  the  gravity  vector  relative  to  the  surface  of  the  shell 
resulting  in  appreciable  components  of  loading  in  the  radial  and  tangential 
directions  (see  equations  8.6.1,  8.6.2;  and  8.6.5)*  These  components  cause 
larger  membrane  loads  which,  in  turn,  make  more  severe  the  bending  action 
at  the  boundary.  It  should  be  borne  in  mind  that  the  asymmetric  shell 
behavior  varies  either  as  sin  0  or  cos  0  (see  section  8.IO). 

There  are  some  additional  interesting  differences  and  similarities 
between  the  symmetric  and  anti- symmetric  behavior  to  be  observed.  The  force- 
results  ,  ^0  ^  again  coincide  for  the  four  sizes.  The  widths 

V  y  ^ 

of  the  edge  zones  as  shown  by  the  curves  of  M  and  and  Q  are 

r  y  r 

essentially  identical.  However,  unlike  the  symmetric  case,  the  maximum 
values  of  the  bending  quantities  increase  significantly  rather  than  decrease 
with  the  size  of  the  shell.  The  maximum  radial  and  tangential  components 
of  displacement  are  the  same  order  of  magnitude  as  the  normal  component. 

Thus,  the  calculation  of  the  resultant  deflection  of  a  point  on  the  shell 

■X' 

must  consider  the  vector  addition  of  u  ,  v  ,  and  w  . 

The  quantities  and  M*  have  not  been  displayed  because  of 
their  small  magnitudes.  For  example,  for  case  4s,  the  maximum  value  of 
Q*  is  -.0058  and  for  ,  the  maximum  value  is  -.004o4.  Both  of  these 
values  are  at  least  an  order  of  magnitude  smaller  than  the  values  for  the 
other  transverse  shear  and  bending  moments,  respectively. 
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Table  8.I3.I 


Case 

\  ^2 
(inches) (inches) 

B.C-  at 

B.C.  at 

Loading 

Is* 

0  230.4 

---- 

Simple  support 

Symmetric- gravity 

2s 

"  460.8 

tt 

tt 

5s 

"  921.6 

— 

tt 

tt 

^4-S 

"  1382.4 

____ 

tl 

tt 

la 

"  230.4 

tl 

Ant i - syra • - gra v i t y 

2a 

"  460.8 

---- 

tl 

tl 

3a 

"  921.6 

---- 

tt 

tt 

4a 

"  1382.4 

It 

tt 

Note : 

Material  properties  are  E  =  10 

7 

^  p  •  s  •  1  •  ^  ~LJ  “  •  3  ^ 

/^  =  .10  Ih.  per 

cubic  inch. 

1 

Focal  length  f  = 

576  inches, 

shell  thickness  h 

=  1  inch. 

*  Calculated  from  shallov  shell  routine  (see  Chapter  7)* 
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8.l4  Group  II  -  Five  Shells  with  an  Outer  Radius  of  460.8  Inches 

The  results  in  this  section  illustrate  the  behavior  of  five  shells 
all  open  at  the  apex  and  with  the  same  outer  radius  of  460.8  inches  but 
terminating  at  various  inner  radii.  All  the  shells  are  simply  supported 
along  both  edges.  Dimensions  and  other  pertinent  data  are  presented  in 
Table  8.14.1.  The  outer  radius  is  equivalent  to  a  slope  of  ^  =  .40  while 
the  inner  radii  correspond  to  =  .l6^  .20,  .2^,  .yo,  and  .35-  Both  the 
symmetric  and  anti- symmetric  behavior  under  gravity  loads  are  presented 
with  the  former  corresponding  to  a  pointing  angle  0°  and  the  latter 

to  =  90 The  curves,  in  addition  to  the  heading  of  three  groups  of 
numbers  (e.g.  8.14.1),  bear  either  the  letter  s,  signifying  symmetric 
behavior  or  a,  signifying  anti- symmetric  behavior.  Results  are  all  in 
normalized  form  (see  page  *43  of  this  report)  . 

The  deflection  results  shown  in  figure  8.14.1  for  the  behavior 
under  a  symmetric  gravity  loading  show  the  same  trend  as  that  exhibited 
by  the  shallow  shells  studied  in  section  7-12. 3i  namely,  that  the  maximum 
deflections  increase  as  the  span  of  the  shell  (we  shall  define  the  span 
to  be  the  difference  between  the  outer  and  inner  radius  of  the  shell)  gets 
smaller.  This  somewhat  unexpected  deformation  behavior  is  due  to  the 
interaction  of  the  membrane  behavior  in  the  interior  of  the  shell  and  the 
bending  behavior  in  the  edge  zones  at  the  boundaries  (see  section  7-12.3) • 

As  can  be  seen  from  Table  8.14.2,  the  deflections  for  the  purely  membrane 
behavior  increase  with  decreasing  span.  In  order  to  correct  the  inherent 
shortcoming  of  the  membrane  solution,  the  asymptotic  solution  can  be 
visualized  as  superimposing  transverse  shear,  moments  and  in-plane  forces  at 
the  boundaries  onto  the  membrane  solution  so  as  to  meet  the  requirement  of 
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zero  w*  displacement  at  the  edges.  These  edge  loads  induce  inextensional 
behavior  (cf.  references  [45]  )  in  the  interior  of  the  shell  and  it 

is  the  inextensional  behavior  which  causes  the  curves  of  w*  (see  figure 
8.14.1)  to  exhibit  humps  near  the  edges. 

* 

The  curves  of  the  bending  moments  M  and  M  shown  in  figures 

r  y 

8.14.3  and  8.14.4  give  a  clear  illustration  of  the  widths  of  the  edge  zone 
at  the  boundaries  of  the  shell.  It  is  interesting  to  note  that  the 
shapes  and  magnitudes  of  each  set  of  the  bending  moment  curves  are  essen¬ 
tially  the  same  in  the  edge  zone  except  for  case  5s  which  does  not  have 
an  edge  zone.  Ihis  latter  case  exhibits  essentially  "plate"  behavior  in 
that  the  deformation  is  predominantly  bending.  The  curve  of  Q*  in 
figure  8.14.5  for  case  5s  also  indicates  "plate"  behavior. 

The  curves  for  the  anti- symmetric  behavior  in  this  section  are 
self-explanatory.  In  all  cases,  the  magnitudes  of  the  various  parameters 
are  smaller  than  their  counterparts  under  symmetric  loading. 
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Table  8.14.1 
Group  II  Examples 


Case 

^1 

(inches) 

^2 

(inches) 

B.C. 

at 

B.C.  at 

Loading 

Is 

184.52 

460.8 

Simple 

support 

Simple  support 

Symmet  r  i  c  -  gravity 

2s 

250.4 

rt 

« 

It 

tt 

5s 

288.0 

tt 

II 

tt 

II 

4s 

545.6 

1! 

II 

tr 

It 

5s 

405.2 

n 

11 

w 

If 

6s 

___ 

If 

— 

tt 

II 

la 

184.52 

>? 

Simple 

support 

It 

Anti-sym. -gravity 

2a 

250.4 

It 

It 

It 

tt 

5a 

288.0 

ff 

11 

tt 

If 

4a 

545.6 

If 

If 

ti 

ft 

5a 

405.2 

ft 

It 

ti 

ft 

Note : 

Material  properties  are 

E  =  lo'^ 

p.s.i.,  V  -  .5, 

-  .10  lb. 

per  cubic  inch. 

Focal  length  f  =  57^  inches^  shell  thickness  h  =  1  inch. 
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Table  8.14.2 

Membrane  Solutions  for  Cases  Is,  2s,  5s  ^s 


Case  Is 


Case  2s 


Case  5s 


Case  4s 


u* 

V* 

N* 

r 

r 

.16 

0 

-  .411 

.466 

.22 

-.oo4o 

-.407 

.495 

.28 

-.0054 

-.402 

.515 

■  54 

-.004l 

-.595 

.535 

.40 

0 

-.586 

.552 

* 

n 

* 

w 

N* 

r 

r 

.20 

0 

-  .415 

.472 

■25 

-.0051 

-.409 

.495 

.50 

-.oo4o 

-.4o4 

.514 

■35 

- .0051 

-.398 

■531 

.40 

0 

-.391 

.548 

* 

u 

* 

V 

N* 

r 

r 

.250 

0 

-.417 

.480 

.288 

-.0019 

-.4l4 

.498 

■  525 

-.0025 

-.409 

■  513 

.563 

-.0019 

-  .405 

.528 

.400 

0 

-■399 

.542 

* 

u 

V* 

N* 

r 

r 

.500 

0 

-.421 

.489 

•325 

- .0009 

-.419 

.502 

■  350 

-.0012 

-.416 

■  513 

.375 

-.0009 

-  .415 

.523 

.400 

0 

-.409 

.534 

64 


o 

iD 

O 


O 

lO 

b 


o 

o 


o 

m 


o 

CVJ 

d 


o 

o 


o 

b 


o 

O 

o 

o 

<N 

fO 

lO 

O 

o 

o 

o 

♦  a> 

z 


65 


0  04  0  08  012  0  16  0  20  0  24  0  23  0.32  0  36  0  40 


M 


0  12  0  16  0  20  0  24  0  28  0  32  0  36  0  40 

r 


66 


5s 

012  0  16  0  20  0  24  0  28  0  32  0  36  0  40 

/ 


67 


3-71  -17881 

FIGURE  8.14.4 
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8.15  Group  III  -  Five  Shells  with  an  Outer  Radius  of  69I.2  Inches 


The  results  in  this  section  illustrate  the  behavior  of  five  shells 
with  the  same  outer  radius  of  6$1.2  inches  but  terminating  at  various  inner 
radii.  Dimensions  and  other  pertinent  data  are  presented  in  Table  8.15.1. 
The  outer  radius  is  equivalent  to  -  .60  while  the  inner  radii  correspond 
to  a  shell  closed  at  the  apex  and  to  slopes  of  =  .20,  .50  and  .55. 

Both  the  symmetric  and  anti- symmetric  behavior  under  gravity  loads  are 
presented  with  the  former  corresponding  to  a  pointing  angle  of  =0®  and 
the  latter  to  -  ^0° .  The  curves^  in  addition  to  the  heading  of  three 
groups  of  numbers  (e.g.  8.15.1)^  bear  either  the  letter  s^  signifying 
symmetric  behavior  or  a,  signifying  anti- symmetric  behavior.  Results  are 
all  in  normalized  form  (see  page  ^^of  this  report). 

In  general,  the  results  are  comparable  to  those  shown  in  section 
8.l4  and  the  same  discussion  is  pertinent.  However,  unlike  the  results  of 
the  previous  section,  the  anti- symmetric  behavior  for  the  larger  span 
shells  (cases  1  and  2)  of  this  group  results  in  larger  values  for  the 
force  resultants,  moments,  and  deflections  than  for  the  symmetric  behavior. 
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Table  8.I5.I 
Group  III  Examples 


3a  se 

(inches) 

R2 

(inches) 

E.C.  at 

B.C.  at 

Loading 

Is 

— 

691.2 

— 

Simple  support 

Sym. -gravity 

2s 

230.4 

tl 

Simple  support 

II 

tt 

5s 

460.8 

!l 

tt 

tt 

4s 

576 

tl 

til 

tl 

tt 

5s 

633.6 

W 

M 

If 

It 

la 

— 

M 

— 

91 

2a 

230.4 

« 

Simple  support 

"  Anti-sym. -gravity 

3a 

460.8 

It 

It 

tt 

tl 

4a 

576 

tl 

91 

tl 

9t 

5a 

633.6 

ft 

t! 

tl 

It 

7 

Note;  Material  properties  are  E  =  10  p.s.i.,  V  -.3)  l^o  ~  *10  lb. 

per  cubic  inch. 


Focal  length  f  »  57^  inches,  shell  thickness  h  =  1  inch. 
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8.l6  Group  IV  -  Seven  Shells  with  an  Outer  Radius  of  921.6  Inches 

The  results  in  this  section  illustrate  the  behavior  of  seven 
shells  with  the  same  outer  radius  of  921.6  inches  but  terminating  at 
various  inner  radii.  Dimensions  and  other  pertinent  data  are  summarized  in 
Table  8.I6.I.  The  outer  radius  is  equivalent  to  "5"  =  .80  while  the  inner 
radii  correspond  to  a  shell  closed  at  the  apex  and  to  slopes  of  ^  =  .52, 
.40,  .50,  .60,  .66,  and  .75.  Both  the  symmetric  and  anti- symmetric  behavior 
under  gravity  loads  are  shown  with  the  former  corresponding  to  a  pointing 
angle  of  =0®  and  the  latter  to  =90®.  The  curves,  in  addition  to 

the  heading  of  three  groups  of  numbers  (e.g.  8.I6.I),  bear  either  the 
letter  s,  signifying  symmetric  behavior,  or  a,  signifying  anti- symmetric 
behavior.  Results  are  all  in  normalized  form  (see  page  4 3 of  this  report). 

In  general,  the  Group  IV  results  are  comparable  to  the  Group  II 

and  III  results  for  shells  with  smaller  outer  radii  and  to  the  Group  V 

results  for  shells  with  a  larger  outer  radii  (read  discussion  in  section 

8.l4).  The  deflections  are  larger  for  the  Group  IV  shells  than  for  the 

smaller  Group  II  and  III  shells.  Note,  however,  that  the  quantities 

associated  with  the  bending  behavior  such  as  edge  zone  width,  M*  ,  M*  , 

X*  y 

,  0^  ,  etc.  have  substantially  the  same  magnitudes. 
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Table  8.16.I 


Group  IV  Examples 


Case  R2  B.C.  at  R^  B.C.  at  R2  Loading 

(inches)  (inches) 


Is  568.64  921.6 

2s  460.8 

5s  576 

4s  691.2  " 

5s  806. 4o  " 

6s  864.00  " 

Ts 

la  568.64 

2a  560.8 

5a  576 

4a  691.2 

5a  806. 4o  " 

6a  864.00  " 

7a  — -  " 


Simple  support  Simple 

tt 


tt 


tt 


tt 


tt 


Simple  support 

tt 


tt 


tt 


tt 


tt 


support  Sym . - gravity 


Anti  Sym. -gravity 


n 

Note:  Material  properties  are  E  =  10  p.s.i.,  ^  =  .3^  ^-.10  lb, 

per  cubic  inch. 

Focal  length  f  »  576  inches,  shell  thickness  h  =  1  inch. 
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8.17  Group  V  -  Ten  Shells  with  an  Outer  Radius  of  l^Q2.k  Inches 

The  results  in  this  section  illustrate  the  "behavior  of  seven 
shells  with  the  same  outer  radius  of  1582.4  inches  but  with  different  inner 
radii.  Dimensions  and  other  pertinent  data  are  summarized  in  Table  8.17.I. 
The  outer  radius  is  equivalent  to  X  =1.20  while  the  inner  radii  correspond 
to  a  shell  closed  at  the  apex  and  to  slopes  of  V  =  .30,  .4o,  .48,  .60,  .75; 
.90,  1.05,  1.075  and  I.15.  Both  the  symmetric  and  anti- symmetric  behavior 
under  gravity  loads  are  shown  with  the  former  corresponding  to  a  pointing 
angle  of  =0°  and  the  latter  to  ^  =  90°.  The  curves,  in  addition  to 
the  heading  of  three  groups  of  numbers  (e.g.  8.17.I),  bear  either  the 
letter  s,  signifying  symmetric  behavior  or  a,  signifying  anti- symmetric 
behavior.  Results  are  all  in  normalized  form  (see  page  4 5 of  this  report). 

In  general,  the  Group  V  results  exhibit  the  same  trends  as  the 
results  for  Groups  II,  III,  and  IV  which  are  shells  with  smaller  outer 
radii.  The  deflections  are  larger  for  the  Group  V  shells  because  the 
outer  radius  is  the  largest  of  all  the  groups.  Note,  however,  that  the 
quantities  associated  with  bending  behavior  such  as  edge  zone  widths  , 

Q*  ,  etc.  do  not  vary  appreciably  between  the  four  groups. 
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Table  8.I7.I 
Group  V  Examples 


Case 

^1 

(inches) 

^2 

(inches) 

B.C.  at 

B.C.  at  Rg 

Loading 

Is 

552.96 

1382.4 

Simple  support 

Simple  support 

Sym. -gravity 

2s 

691.2 

U 

tl 

It 

ft 

3s 

864 

It 

It 

ri 

tl 

4s 

1036.8 

It 

M 

It 

tl 

5s 

1209.6 

It 

It 

It 

tl 

6s 

1296.0 

tl 

tl 

It 

tl 

7s 

— 

II 

— 

It 

tl 

la 

552.96 

It 

Simple  support 

II 

Anti-Sym. -gravity 

2a 

691.2 

It 

It 

ri 

II 

3a 

864 

M 

It 

It 

II 

4a 

1036.8 

It 

It 

It 

II 

5a 

1209.6 

It 

It 

It 

II 

6a 

1296.0 

tl 

It 

II 

II 

7a 

— 

tl 

— 

It 

It 

8a 

345.6 

It 

Simple  support 

It 

tl 

9a 

460.8 

It 

It 

M 

tl 

10a 

1324.8 

It 

II 

tl 

tl 

Note;  Material  properties  are  E  -  10  p.s.i.,  V  -  .J),  -  .10  lb. 

per  cubic  inch. 

Focal  length  f  =  576  inches,  shell  thickness  h  =  1  inch. 
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8.l8  Group  VI  -  Four  Shells  Which  are  Free  at  the  Inner  Boundary 

The  results  in  this  section  illustrate  the  behavior  of  four 
shells  with  the  same  outer  radius  of  i^^0.8  inches  but  with  different  inner 
radii.  In  this  section,  each  shell  is  simply  supported  along  the  outer 
edge  and  is  free  at  the  inner  edge.  Dimensions  and  other  pertinent  data 
are  summarized  in  table  8.I8.I.  The  outer  radius  is  equivalent  to  =  .4 
while  the  inner  radii  correspond  to  slopes  of  y  =  .16,  .20,  .25,  and  .50. 

w  •X’ 

Only  w  and  have  been  shown. 

It  is  observed  that  the  deflections  for  the  symmetric  loading 
increase  as  the  span  of  the  shell  decreases.  This  is  contrary  to  the  usual 
behavior  of  planar  structures  and  dramatically  illustrates  the  manner  in 
which  a  shell  derives  its  rigidity.  The  missing  part  of  the  shell  (i.e., 
the  part  which  was  removed  to  make  the  hole)  supplies  support  to  the 
remainder  of  the  shell.  (One  can  visualize  what  would  happen  to  a  simple 
cable  if  a  piece  of  it  were  removed.)  Since  the  shell  is  a  two  dimensional 
structure,  it  will  not  collapse  because  the  loss  of  the  cable  action”  is 
compensated  by  an  increase  in  the  hoop  action.  This  mode  of  behavior  is 
not  as  efficient  as  the  ”cable  action”  and  hence,  leads  to  larger  dis¬ 
placements.  The  trend  of  increasing  deflection  with  decreasing  span  does 
eventually  reverse  (as  it  must)  but  at  a  span  which  corresponds  to  a  very 
narrow  ring  (see  shallow  shell  results).  Note  that  the  bending  moment 
M*  exhibits  a  rather  slow  variation  with  changing  span. 

In  general,  the  results  are  the  same  as  those  previously  shown 
for  shallow  shells  (see  sections  7*12.4  and  7*12.5)* 
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The  anti- symmetric  case  does  not  exhibit  the  same  trend  of  in¬ 
creasing  deflections  with  decreasing  span.  In  the  anti- symmetric  case, 
the  larger  components  of  the  gravity  load  vector  are  aligned  along  the 
tangent  plane  rather  than  along  the  surface  normal.  This  means  that  there 
is  less  bending  action  in  the  anti- symmetric  loading,  and  it  is  the 
bending  action  which  leads  to  the  results  which  have  been  discussed. 
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Table  8.I8.I 
Group  VI  Examples 


Loading 


R^  B.C.  at  R^ 


(inches) 

(inches) 

Is 
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Note:  Material  properties  are  E 
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8.19  Group  VII  -  Five  Shells  which  are  Free  at  the  Inner  Boundary 

The  results  in  this  section  illustrate  the  "behavior  of  four 
shells  with  the  same  outer  radius  of  921.6  inches  hut  with  different  inner 
radii.  As  in  section  8.I8,  each  shell  is  simply  supported  at  the  outer 
radius  and  is  free  at  the  inner  radius.  Dimensions  and  other  pertinent 
data  are  summarized  in  table  8.19.1.  The  outer  radius  is  equivalent  to 
V  =  .8  while  the  inner  radii  correspond  to  slopes  of  Y  =  .32,  .40,  .50^ 
.60,  and  .70.  Only  w  and  have  been  shown.  Results  are  presented  in 
normalized  form. 

The  Group  VII  results  are  similar  to  those  of  Group  VI  and  the 
discussion  in  section  8.I8  is  applicable. 
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Table  8.19.I 
Group  VII  Examples 


Case 

^1 

( inches) 

Is 
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806.4 
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Rg  B.C.  at 

(inches) 

921.6  free 


B . C .  at  R2 
Simple  support 


tt 
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Note:  Material  properties  are  E*10p.s.i.,  J=.5 
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8.20,  6.21,  6.22  Group  VIII  -  Asymmetric  Behavior  of  Three  Shells  Simply 

Supported  at  the  Outer  Boundary 


In  previous  sections,  the  symmetric  behavior  and  anti- symmetric 
behavior  which  correspond  to  pointing  angles  of  0°  and  90°  were  pre¬ 
sented.  In  this  section,  the  deflection  patterns  for  the  shell  at  pointing 
angles  of  ^  =  15°,  30° ^  ^5°^  60°,  and  90°  are  shown  in  the  form  of  polar 

plots . 

As  has  already  been  discussed  (see  section  7*12.7  of  part  III) 
large  antennae  will  generally  be  erected  in  the  face-up  position  (  ^  0°) 

Thus,  it  becomes  of  interest  to  determine  the  deflections  of  the  antenna 
relative  to  the  face-up  position.  The  displacements  normal  to  the  middle 
surface  for  gravity  loads  are  given  by 


~  (O  cos  ip  +  w  (  ^  )  sin  ^  sin  0 

S  Q> 


vhere 

w.  ( n 


8.20.1 

symmetric  part  of  the  solution 

anti- symmetric  part  of  solution  at  0  =  IL- 

2 


The  plots  in  the  preceding  sections  are  of  and  w  . 

Let  w  (  0  ,  0)  be  the  deflection  relative  to  the  antenna  in  its 
face-up  position,  i.e.,  relative  to  ~  0.  Then 


w  (5^,0)  =  w(5'^0)-w  8.20.2 

s 

which  leads  to 

w  (  K ,  9)  =  w^(cos  ^  -  1)  +  w^  sin  sin  0  8  20  3 

The  three  shells  analyzed  in  this  section  have  the  same  dimen¬ 
sions  as  shells  previously  considered.  Dimensions  and  other 
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pertinent  data  are  as  follows: 

7 

E  -  10'  p.s.i. 

V-.  .3 

.10  lb.  per  cubic  inch 
focal  length  f  ^  5T6  inches 

shell  thickness  h  =  1  inch 

All  three  shells  are  closed  at  the  apex  and  the  boiindaries  correspond  to 
radii  of  460.8  inches,  691-2  inches  and  921-6  inches,  respectively. 

The  results  are  presented  in  the  form  of  contour  plots  wherein 
lines  of  constant  "w*  are  shown.  A  polar  representation  has  been  used  in 
which  the  horizontal  radius  is  the  coordinate  on  the  shell  given  by  0  =  0. 
Radii  emanating  from  the  origin  are  lines  of  negative  0  below  the  horizon¬ 
tal  and  positive  0  above  the  horizontal.  Each  plot  constitutes  the 
deflected  shape  of  the  shell  in  the  normal  direction  relative  to  the 
face-up  position,  i.e.,  of  "w*  ,  for  various  pointing  angles  ^  .  Since 

the  behavior  is  symmetric  about  a  vertical  plane,  only  the  region  of  the 
shell,  ^0  ^  ^  has  been  shown. 

The  curves  are  self-explanatory.  Each  contour  is  roughly  in 
the  shape  of  a  U  with  the  legs  approximately  horizontal.  Positive  values 
mean  that  the  shell  has  deflected  inward.  Note  that  the  largest  values 
are  positive  and  occur  in  the  upper  half  of  the  shell. 
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FIGURES  8.20.4,  5,  &  6 
ASYMMETRIC -VARYING  i|i 


139 


140 


141 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

(VJ 

ro 

in 

U) 

o 

iO 

in 

ro 

<VJ 

d 

d 

d 

d 

d 

d 

ro 

d 

d 

d 

d 

d 

o 

d 

o  » 

1 

1 

1 

1 

1 

II 

o 

O 

fO 

w 

d 

6 

o 

o 

o 

o 

o 

<M 

o 

o 

U) 

o 

CO 

00 

(D 

(Vi 

d 

d 

d 

d 

d 

d 

d 

d 

o 

1 

1 

1 

1 

142 


ASYMMETRIC  -  VARYING  if. 
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